Micro-magnets for coherent control of spin-charge qubit in lateral quantum dots 
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A lateral quantum dot design for coherent electrical manipulation of a two-level spin-charge system 
is presented. Two micron-size permanent magnets integrated to high-frequency electrodes produce 
a static slanting magnetic held suitable for voltage controlled single qubit gate operations. Stray 
held deviation from the slanting form is taken into account in the Hamiltonian describing the two- 
level system, which involves hybridization of a single electron spin to the quantum dot's orbitals. 
Operation speed and gate fidelity are related to device parameters. Sub 100 ns ir pulse duration 
can be achieved with lattice fluctuations coherence time of 4 ms for GaAs. 
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Quantum dots (QDs) are man-made structures which 
can confine conduction electron in semiconductor to a 
nanometer size volume [U, H[ . ln the few-electron regime, 
the spin of individual electrons, which is a natural two- 
level system, can be used as a quantum bit for imple- 
menting scalable quantum computing Q . In this context 
various experiments have been performed on robustness 
of electron spin in single- and double-dot systems 0, Q , 
spin correlation [(| 0, H| , electrical readout of single elec- 
tron spin [9fl a nd coherent manipulation of two-electron 
spin states [10[ ■ The rotation of a single electron spin or 
qubit operation, which is the most fundamental quantum 
operation, was recently demonstrated by generatin g an 
on-chip electron spin resonance (ESR) magnetic field [ll| . 
An alternative approach is to modulate a QD electric field 
in a non- uniform magnetic field [12J . This scheme elimi- 
nates the need for an externally applied ac magnetic field, 
and does not require spin-orbit coupling, as opposed to 
earlier work on electron spin control based on g-tensor 
modulation, and on electric fields [HI, [l4| ■ Instead, ESR 
is achieved by applying microwave gate voltage pulses, 
letting the electron position oscillate in a static slanting 
Zeeman field. This effectively provides the electron spin 
with the necessary time-dependent transverse magnetic 
field. 

We present in this paper a lateral quantum dot design 
to perform electrical ESR in a slanting Zeeman field. The 
design involves a realistic permanent magnets configura- 
tion to produce the required inhomogeneous static mag- 
netic field and high-frequency electrodes for fast charge 
control. The geometry used is that of well-established 
few-electron quantum dots defined in semiconductor two- 
dimensional electron gas (2DEG) by surface gates [l5[. 
The coupling between the spin and orbital degrees of free- 
dom allows fast qubit rotations with decoherence time 
only slightly reduced from the intrinsic spin dephasing 
time for GaAs quantum dots. 

The device layout is presented in Fig. [T]a. Four metal- 
lic gates, located on the surface of a semiconductor het- 
erostructure, form Schotky contacts to a 2DEG located 
at a distance d below. For GaAs/AlGaAs wafer, d is 
typically 100 nm. By applying suitable dc gate bias, 



a quantum dot confining a single excess electron, sep- 
arated from source and drain reservoirs by tunnel bar- 
riers, is formed in the 2DEG. The QD confining poten- 
tial is composed of a transverse part V (z) coming from 
the band structure and in-plane part V (x, y) electro- 
statically defined by gates L, P, R and T. For typical 
gate bias, the confinement can be modeled as V(x, y) = 

\rruv 2 x 2 + y 2 + "f(y — jj) 2 [lj| where 7 characterizes 

nonparabolicity and the length D is the diameter of semi- 
circular confinement by the T gate. The slanting mag- 
netic field is produced by two ferromagnetic strips po- 
sitioned on each side of the gate structure (see below). 
Each strip is connected to a high-frequency voltage port 
through impedance-matched on-chip waveguide to enable 
GHz electric field modulation. Anti-phase modulation 
of the strips potential at frequency v produces a uni- 
form electric field E = Eq sin (2-Kvt) x inside the quan- 
tum dot with amplitude Eq proportional to the potential 
drop Vq between the strips and inversely proportional 
to the strip separation s. From numerical simulations 
of the microwave field produced by the strips, we find 
E a = 7.71 x 10~ 7 V /s in S.I units. 

The external magnetic field, required for the ESR 
scheme, is applied in the plane of the QD whose strength 
Bo (few Teslas) is much higher than the coercive field of 
the patterned magnets (ranging from zero to hundreds 
of Gauss). In these conditions, the strips are magnetized 
uniformly in the direction parallel to the external field. 
The magnetic moment per unit volume is equal to the sat- 
uration magnetization M of the ferromagnetic material. 
The resulting stray magnetic field is illustrated in Fig. 
[T]a. In the magnets plane (z = d + t/2), the stray field is 
parallel to the external field. In the 2DEG plane (z = 0) 
however, a perpendicular component is present due to the 
curvature of the field lines. Figure [TJb shows the profiles 
of the parallel (B^. 1 ) and transverse (B^) components 
of the stray field at the 2DEG location. The transverse 
component peaks with opposite sign at the strips edges 
x = ig/2 while the in-plane component switches sign at 
these positions. These properties produce, near the QD 
location, a total magnetic field of the form 
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B=[B Q + SB + b SL z + a (z 2 



3bzx 2 — fez 3 ] x + \psLX + 2axz — 3bxz 2 + fex 3 ] I 
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which satisfies Maxwell's equations V-B = and VxB = 
0. 

The terms in Eq. ([1]) up to the first order in x and 
z displacements correspond to the slanting form previ- 
ously considered in Ref. [12]. These terms are described 
by a uniform shift SBq of the external field and gra- 
dients with slope bsL on the order of lT/fim. Coef- 
ficients a and fe give respectively the second and third 
order corrections to the slanting form. For displace- 
ment smaller than ±25 nm, we find that Eq. ([T]) fits 
very well the magnetic field profile obtained by numer- 
ical simulation [2l| . with departure of less than 0.1%. 
Because of the two-dimensional nature of the QD (i.e. 
strong transverse confinement), the effect of the z terms 
in Eq. {J) on the electron dynamics can be averaged 
out. This procedure involves averages (z fc ) = (£o| z k |£o) 
over the ground state wavefunction £o (z) of V (z). Term 
k = 1 is equal to zero and terms k = 2, 3 renormalize 
the shift and slope to SB = SBq + a(z 2 ) — fe(z 3 ) and 

1>sl = bsL — 3fe(z 2 ) to give the effective QD magnetic 
field (b) = 



B + SB - ax 2 



b.SLX + bx 3 



We now relate the proposed device to the control of 
the hybrid spin-charge qubit. With M = 0, the electron 
spin confined into the QD is described by the Hamilto- 
nian H = E + (pi + P 2 y ) fern + V (x, y) - gfi B B a x /2 
where Eq is the ground state energy of the transverse 
Hamiltonian, m the effective mass, g the g-factor and [is 
the Bohr magneton. Throughout the text, we use mate- 
rial parameters of GaAs: m = 0.067m e and g = —0.44 
where m e is the electron mass. The Pauli spin ma- 
trices a = {cr x ,a y ,a z ) are decoupled from the QD or- 
bitals <p n (x,y). The eigenenergies and eigenfunctions of 
Ho are e na = Eq + e n — ^g/j,BB a with eigenfunctions 
(x,y,z\n,a) = £ (z) 4> n (x, y) ij} a where n = 0,1,2,..., 
a = ±1 and ip a is the spinor of a x eigenstates. The 
ground state wave function is spin split by the Zeeman 
energy which is assumed to be smaller than orbital exci- 
tation energy: g^tsBo < ei — eo. 

With M > 0, the Hamiltonian becomes H — Hq + 
W x (x) a x + W z (x) o z with position dependent pertur- 
bations W x (x) = —gfJ-B [SB — ax 2 ) /2 and W z (x) = 



~9P"B (J>slx + bx 3 ^j /2. The relatively large value of 5Bq 

shifts the Zeeman energy of the QD spin from the nearby 
reservoirs, a useful feature for read-out schemes based 
on spin to charge conversion [13 . Terms W x , z (x) mix 
the electron spin states \a) with the orbital states \n). 
For example, the two lowest energy states, which con- 
stitute our "hybridized" spin-charge qubit, are, for the 
harmonic potential (7 = 0), \ga) w Co a \n x — 0, a) + 
Ci- a \n x = 1, — cr) where n x =0,1 corresponds to the 



first two orbitals of the one dimensional harmonic po- 
tential (of the x coordinate) [12]. For the magnetic field 
profile shown in Fig. l.b and typical confinement energy 
hujo = lmeV, we find Ci-cr/Coa ~ 0.0003. The mixing 
between the spin and orbital states is therefore relatively 
small. Neglecting for the moment possible misalignment 
between the magnetic strips and the quantum dot, only 
the off-diagonal matrix elements of the high-frequency 
QD electric potential operator —eEgx sin(2nvt) remain 
in the subspace spanned by {\ga — +1) , \gu = — 1)}. 
The Hamiltonian of the qubit is thus of the ESR form 
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where h is the Planck constant and a x and a x are the 
Pauli matrices of the effective two-level system. 

The Larmor frequency vqi (hvoi corresponding to the 
qubit energy level separation) is slightly smaller than 
unperturbed Zeeman energy g^isiBo + SBq). The ESR 
transverse field e x , which controls qubit rotation speed, 
is proportional to the modulation voltage Vq and magne- 
tization M. By adjusting the phase of the modulation, 
arbitrary qubit rotations can be achieved when v = vq\. 
The predicted dependence of Pqi and n pulse duration 
time, tv = h/e x , on strips's thickness t, separation s 
and Vq is presented in Fig. O Results are obtained by 
least square fit of the simulated stray magnetic field and 
microwave electric field followed by exact diagonaliza- 
tion of the corresponding spin-charge Hamiltonian H for 
Bq = 2 T. In the calculation, we use confining potential 
parameters Hloq = lmeV, 7 = 1, and D = 2.9y // h/mujQ 
and magnetization of cobalt hqM = 1.8 T. For a given 
strip separation, the thicker are the strips, the bigger are 
^01 and e x . For thickness greater than approximately 
4d, Uqi and tv start to saturate to their t — > 00 value. 
Only small improvement is therefore achieved in using 
strips thicker than 0.4 fim. Small s should be used as 
it improves drastically qubit rotation speed. For rela- 
tively small electric field Eo w 5mV//mi (correspond- 
ing to point Vo = 3mV in inset of Fig. 2a), the ro- 
tation speed is comparable to the fastest Rabi oscilla- 
tions achieved recently in GaAs lateral quantum dot with 
on-chip ESR magnetic coil operating at much lower fre- 
quency (v = 200MHz) [ll|- We speculate that for the 
GHz regime, Rabi oscillations faster than with conven- 
tional ESR should therefore be observed. The misalign- 
ment between the strips and the QD introduces an oscil- 
lating term e z sin (2tiv£) a z in the qubit Hamiltonian. For 
misalignment as large as 100 nm, we find that this term 
does not influence gate operations, i.e. the gate fidelity 
for 7T pulse remains one with error less than 10 . 

Usual spin qubit, in the absence of spin-orbit coupling, 
are naturally protected against lattice fluctuations. For 
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our case, because the spin states are coupled to the elec- 
tron's orbitals, effect of phonons has to be considered. As 
done by Tokura et al. [l2| for one-dimensional systems, 
acoustic phonon scattering causes orbital relaxation be- 
tween the qubit states. Slightly modifying the previous 
analysis to our two-dimensional quantum dot, we esti- 
mate a relaxation time Xi ~ 2.1 msec at Bq = 2 T, which 
is dominated by transverse piezoelectric scattering. Since 
the dephasing effect without relaxation is negligible, the 
coherence time is 27\. The corresponding quality factor 
Q is estimated by coherence time divided by tv , which is 
order of 10 4 , similar to the one-dimensional system con- 
sidered previously. Q factor of the same order has been 
estimated for the case of spin-orbit coupling [HI]. The 
main source of decoherence comes from the hyperfine in- 



teraction with GaAs nuclei spins. The nuclei will induce 
fast dephasing (with time averaged coherence time of 10- 
20 ns) as demonstratred by the work of Petta et al [l(| • 
However, because the hyperfine field fluctuates on a slow 
time scale compared to gate operation speed, this effect 
can be corrected using compensation [Hj, [HI or projec- 
tion techniques [2(3]. 

We acknowledge W.G. van der Wiel for useful dis- 
cussions and T. Taniyama and Y. Sekine for inputs on 
micro-magnets fabrication. ST thanks financial supports 
from the DARPA (DAAD19-01-1-0659) of the QuIST 
program, the Grant-in- Aid for Scientific Research A (No. 
40302799), SORST-JST and IT Program, MEXT. 



[1] L. P. Kouwenhoven and C. M. Marcus, Phys. World 
(June 1998). 

[2] L. Jacak, P. Hawrylak, and A. Wojs, Quantum Dots 

( Springer- Verlag, Berlin, 1998). 
[3] D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 

(1998). 

[4] T. Fujisawa, D. Austing, Y. Tokura, Y. Hirayama, and 

S. Tarucha, Nature 419, 278 (2002). 
[5] J. R. Petta, A. C. Johnson, A. Yacoby, C. M. Marcus, 

M. P. Hanson, and A. C. Gossard, Phys. Rev. B 72, 

161301 (2005). 

[6] M. Ciorga, A. Wensauer, M. Pioro-Ladriere, M. Ko- 
rkusinski, J. Kyriakidis, A. S. Sachrajda, and P. Hawry- 
lak, Phys. Rev. Lett. 88, 256804 (2002). 

[7] S. Sasaki, S. Amaha, N. Asakawa, M. Eto, and 
S. Tarucha, Phys. Rev. Lett. 93, 017205 (2004). 

[8] T. Hatano, M. Stopa, and S. Tarucha, Science 309, 205 
(2005). 

[9] J. M. Elzerman, R. Hanson, L. H. W. V. Beveren, 
B. Witkamp, L. M. K. Vandersypen, and L. P. Kouwen- 
hoven, Nature 430, 431 (2004). 

[10] J. Petta, A. Johnson, J. Taylor, E. Laird, A. Yacoby, 
M. Lukin, C. Marcus, M. P. Hanson, and A. Gossard, 
Science 309, 2180 (2005). 

[11] F. H. L. Koppens, C. Buizert, K. J. Tielrooij, I. T. Vink, 
K. C. Nowack, T. Meunier, L. P. Kouwenhoven, and 
L. M. K. Vandersypen, Nature 442, 766 (2006). 

[12] Y. Tokura, W. G. van der Wiel, T. Obata, and 
S. Tarucha, Phys. Rev. Lett. 96, 047202 (2006). 

[13] Y. Kato, R. C. Myers, D. C. Driscoll, A. C. Gossard, 
J. Levy, and D. D. Awschalom, Science 299, 1201 (2003). 

[14] E. I. Rashba and A. L. Efros, Phys. Rev. Lett. 91, 126405 
(2003). 

[15] M. Ciorga, A. S. Sachrajda, P. Hawrylak, C. Gould, 
P. Zawadzki, S. Jullian, Y. Feng, and Z. Wasilewski, 
Phys. Rev. B 61, R16315 (2000). 

[16] J. Kyriakidis, M. Pioro-Ladriere, M. Ciorga, A. S. 
Sachrajda, and P. Hawrylak, Physical Review B (Con- 
densed Matter and Materials Physics) 66, 035320 (2002). 

[17] H.-A. Engel and D. Loss, Phys. Rev. Lett. 86, 4648 
(2001). 

[18] A. V. Khaetskii and Y. V. Nazarov, Phys. Rev. B 64, 
125316 (2001). 



[19] J. M. Taylor, H.-A. Engel, W. Ur, A. Yacoby, C. M. 
Marcus, P. Zoller, and M. Lukin, Nature Physics p. 177 
(2005). 

[20] D. Klauser, W. A. Coish, and D. Loss, Physical Review 

B 73, 205302 (2006). 
[21] We used Mathematica Radia package available at 

http://www.esrf.fr/ for stray field calculations. 



4 




-0.5 0.5 

x(|im) 



FIG. 1: (a) Top and cross-sectional views of the lateral quan- 
tum dot device consisting of four metallic gates (yellow) and 
two ferromagnetic strips (blue) patterned at the surface of a 
semiconductor heterostructure. The thick arrows indicate the 
direction of the external magnetic field Bo and magnetization 
M. The origin is fixed to the quantum dot position, (b) In- 
plane (-B* 7 ) and transverse (B^ 1 ) profile of the stray magnetic 
field produced for two Co micromagnets (/ioM = 1.8 T) cal- 
culated at z — 0. Dimensions d, t, s, w and I are set to 0.1, 
0.2, 0.5, 0.5 and 2.0 fim. 
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FIG. 2: Calculated dependence of (a) spin-charge qubit's Lar- 
mor frequency voi and (b) gate operation speed tv for ESR 
voltage Vo = 1 mV on magnetic strips's thickness t and sepa- 
ration s. Black, red, green, blue, yellow and magenta curves 
correspond to value s of 0.5, 0.6, 0.7, 0.8, 0.9 and 1.0 /im 
respectively. Other dimensions are same as in Fig. [1] Mag- 
netization noM = 1.8 T is assumed. Inset: Rabi frequency 
frabi = 1/(2t t ) as a function of Vo for t = 0.2 and 
s = 0.5 /jm. 



